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Gn : Abstract 

Dilepton spectra are calculated within the microscopic transport model 
UrQMD and compared to data from the CERES experiment. The invariant 
mass spectra in the mass region between 300 and 600 MeV depend strongly 
on the mass dependence of the p meson decay width which is not sufficiently 
determined by the Vector Meson Dominance (VMD) model. A consistent ex- 
planation of both the recent Pb+Au data and the proton induced data can 
be given without additional medium effects. 
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I. INTRODUCTION 



Recently, electromagnetic radiation in form of lepton pairs has been observed at CERN 
in the CERES experiment. A strong enhancement above the cocktail of hadronic decays has 
been reported. This is interesting, because dileptons may escape nearly undisturbed from 
the hot and dense zone which is formed in the heavy ion collision. Dileptons could probe 
the intermediate stage of these interactions. On the other hand, hadrons are thought to 
probe the later freeze-out stages. This folklore may not be entirely correct: A number of 
investigations indicate that hadronic signals such as strangeness and entropy, as well as flow 
are developped in the early phase, while a large fraction of the observed dileptons is simply 
due to post-freeze-out decays ("feeding") of isolated hadrons. 

The fact that also dileptons are emitted from various sources and during the entire reac- 
tion makes the interpretation of this signal not as easy as hoped for. Electromagnetic probes 
are definitely not a pure signal from the most interesting high density stage only! Rather, 
their spectrum is composed of different contributions integrated over the full interaction 
dynamics. Transport models need to be applied to decompose the dilepton cocktail into its 
various individual contributions. 

The most discussed feature of the data is the enhancement of intermediate mass lepton 
pairs (300 MeV < M ee < 700 MeV) as compared to the hadronic cocktails based on primary 
collisions. While the enhancement at BEVALAC energies is still not understood Hp-fl; first 
predictions for Pb+Au collisions from transport models including secondary production of 
dileptons (e.g. — > e + e~) indicated only a small (if any) deviation from the CERES data 
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|5|||. A large variety of attempts has been made to explain this deviation. The two most 
popular of them are the hypothesis of lowering vector meson masses (the Brown- Rho scaling) 
and the collisional broadening of the spectral functions in the hadronic environment 

mm- 

Here we analyse the recent CERES data on the multiplicity and pt of the lepton pairs. 
First, we give a survey of the dilepton production in the framework of the ultrarelativistic 
quantum molecular dynamics model, UrQMD |Tj|. The present paper demonstrates that a 
state-of-the-art transport model based on hadronic and string degrees of freedom can predict 
the present data without assuming phase transitions or unconventional in-medium effects. 
This compares well to recent results from a simple transport model ||13| . 



II. DESCRIPTION OF THE MODEL 

Dilepton production in the framework of the UrQMD model has been considered in some 
detail in ref. Q. In this model dileptons are produced perturbatively, mainly in hadronic 
decays. Dalitz decays of neutral mesons and direct decays of vector mesons are taken into 
account. The incoherently summed pn bremsstrahlung as well as the A(1232) Dalitz decay 
are of minor importance at CERES energies because the system is dominated by mesonic 
degrees of freedom. All dilepton decay widths have been specified in H. 

For the vector meson decay width Y p _> e + e - (M) we use the typical M -3 dependence which 
results from vector meson dominance (VMD model). There seems to be some inconsistency 
between the original formulation of VMD (see e.g. Sakurai flip, which results in a M~ 3 
dependence for T p ^ e + e - and the extended VMD from Kroll, Lee and Zumino [15] implying a 



dilepton decay width proportional to M (see e.g. ref. |1|). The discrepancy between these 
two approaches seems to be as M 4 ! However, if one considers the channel hit — > e + e~ , the 
extended VMD requires also a direct coupling between the electromagnetic and the hadronic 
channel. If this term is included (as well as the interference), one ends up with the overall 
M~ 3 dependence again. 

In the extended VMD the direct term can be motivated only by a TTTc-f coupling. There- 
fore it is not clear how to treat the dilepton decay of p mesons which themselves come from 
heavy resonance decays or even string fragmentation. Similar problems arise for the u and 
4> decays. However, due to their small total widths, this effect does not manifest itself in the 
dilepton mass spectra and is thus not discussed separately here. In the UrQMD model the 
following procedure is applied: 

For p mesons which are produced by 7T7r-annihilations or in decays of baryonic or mesonic 
resonances the dilepton decay width scales like M -3 . The dilepton width for p mesons from 
string fragmentation is chosen according to the direct term of the extended VMD model 
which scales like M. Other approaches are possible and yield different results: see discussion 
below. 



III. RESULTS 

The undisturbed dilepton mass spectrum for proton on beryllium collisions at 450 GeV 
incident energy is shown in fig. |l[ This 'naked' spectrum could be verified only with an ideal 
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detector covering the Air phase space angle and having a perfect momentum resolution (the 
bin width is 8 MeV). It shows all details of the model calculation. The CERES data |l7j are 
drawn for orientation only and should not be directly compared to the actual calculation. 
One sees, however, that roughly two orders of magnitude are lost due to the acceptance. 

At low invariant masses the overwhelming background is due to the Dalitz decays of 
the long-living mesons tt° and rj. Between 500 and 750 MeV the u Dalitz and the direct 
p meson decays dominate the spectrum. At higher energies the characteristic peaks of the 
direct uj and <fi decays overlap the p decay. The rj' Dalitz decay or baryonic channels like 
pn bremsstrahlung play only a subordinate role. The asymmetrical shape of the p mass 
distribution is due to phase space limitations and especially due to the M~ 3 scaling of the 
decay width. At invariant masses around 1.2 GeV there is a drop in the p distribution. This 
is due to a limitation of the model not to create resonances from string fragmentation with 
masses exceeding the pole mass more than three widths (M < Mr + 3r#). 

For comparison with the experimental data from the CERES collaboration, one needs to 
correct the model calculations for the limited detector acceptance and momentum resolution. 
Because the detector covers only a small geometric region, every detected particle must have 
a certain angle to the beam axis. This requires the pseudo rapidity of the lepton r] e ± to be 
2.1 < rj e ± < 2.65. A cut on small opening angles O e + e - > 0.035 mrad of the electron-positron 
pair has to be matched because close pairs cannot be identified by the RICH detectors. The 
individual leptons which have passed the initial RICH detectors are then deflected in a 
magnetic field with a momentum resolution parameterised as 

Ap 



V" 2 + m 2 , (i) 
p v 

where a and j3 depend on the detector setup and are listed in table 1 ||18|| . The momenta of 
the electron and positron are washed out accordingly with a Gaussian distribution function 
of width — . Finally, these tracks have to survive a transverse momentum cut, which 
is typically pt > 50 MeV for the proton induced and pt > 200 MeV for the ion induced 
reactions. 

In fig. we present the dilepton mass spectra for proton on beryllium collisions including 
the CERES acceptance and resolution. As compared to fig. [I] the momentum resolution has 
washed out the peak structures and the low mass pairs have been suppressed due to the 
acceptance. As an effect of the finite resolution one can see also contributions from the direct 
p° decay below the two pion threshold. The transport calculation slightly overestimates the 
data around the vector meson poles but seems to be consistent with the overall data set. At 
invariant masses higher than 1.3 GeV the data possibly indicate a need for additional sources 
like direct production in meson+meson collisions (see e.g. |19|) which were not considered 
in the present calculation. 

The UrQMD result for the heaviest system measured by CERES can be found in fig. |3| 
together with the data from the '96 run JIB]. These default calculations do not show the 



typical trend to underestimate the data in the region around 400 MeV. Like in a number of 
other free or conventional transport calculations the data point at 780 MeV, just at the p/u 
mass, is slightly overestimated [|10],|rT],|13j . This has been interpreted that 'some of the signal 
should be distributed away' by medium effects. However, in view of the good agreement 
of the conventional calculation at lower masses one must ask, where this strength should 
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go. Furthermore, we have observed that introducing collision widths results in a significant 
enhancement of vector meson production, which is not compatible to the data. 

Two data points are missed at 140 MeV and 180 MeV. By lowering the p t cut to less 
than 200 MeV this dip can be artificially removed on cost of overshooting the data points 
at 63 MeV and 88 MeV. Our interpretation of this observation is that the UrQMD p t - 
distribution of pions and etas is too soft to be compatible to this data. Furthermore the 
spectral density of the p meson in nuclear matter might extend to below 2m 7r because the 
p may decay into a nucleon nucleon-hole pair. This contribution is neglected in all on-shell 
transport models like UrQMD, but might fill in this gap. 

The results of the transport simulation depend on the details of the VMD model (fig. |]). 
In the default calculation, the dilepton decay width of p mesons scales like M~ 3 if they 
are from hadronic sources and like M for p's from strings (see the discussion above). If no 
VMD was considered at all, i.e. assuming constant dilepton decay widths for the p, a clear 
enhancement of the data over the model calculation can be observed (dashed line in fig. |]). 
However, as pointed out above, neglecting the M~ 3 -dependence of the p meson dilepton 
decay width is not a valid assumption. In the hydrodynamical calculations of ref. |2IJ the 
M~ 3 factor is missed, which might be the reason for their underestimation of the CERES 
data. 

On the other hand, one can assume 'total VDM' by requiring that the vector mesons 
which were produced in string fragmentation scale like the 'hadronic' p's with M~ 3 . This 
leads to an overestimation of the Pb+Au data (see dotted line in fig. [|) and also the p+Be 
data are overestimated. From a theoretical point of view this approach has to be considered 
as an upper limit to the contribution of dilepton emission from p mesons from string frag- 
mentation. Contrary, the default calculation provides a lower limit for a realistic dilepton 
spectrum. The range between both approaches is the result of our uncertainty how to treat 
dilepton production from string fragmentation. 

The recent lead data have also been analysed for the transverse momenta of the lepton 
pairs (i.e. the virtual photons). The idea was to cut on pairs emitted from low-pt vector 
mesons for which the interaction with the surrounding medium should be strongest. Con- 
sequently, in-medium effects would result in a stronger enhancement of the low-p t data set 
over conventional models. The results of the corresponding UrQMD analysis are shown in 
fig. [| in comparison to the CERES data [0 . 

Note, first of all that the dip below M e + e - = 200 MeV which was already observed in 
fig. ^| can be found only in the high p t region (lower frame). This indicates once more that 
UrQMD predicts less high-p 4 7r° and rj mesons. 

In the p region a reasonable agreement between simulation and data can be found both 
for the low-pt (upper frame) and the high-p t sample. In particular the low-pt data are 
consistent with the present UrQMD calculation - without a need for additional medium 
effects. 

IV. CONCLUSION 

In summary, the p+Be and the Pb+Au CERES data can be consistently explained by 
the UrQMD transport simulations. A reasonably good agreement was found especially in the 
region of 0.3 GeV < M < 0.6 GeV. It was discussed that the parameterisation of the vector 
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meson decay width is not constrained from first principles. The variation between different 
approaches for T p ^ e + e -(M) manifests itself in dilepton spectra clearly underestimating or 
slightly overshooting the CERES data. Further efforts in theory are required to provide 
a microscopic scheme for F p _ te + e -(M). Measurements with higher precision could help to 
identify the u> peak in the invariant mass spectra and thus pin down the shape of the p 
meson. In addition more experimental information on the elementary channels p+p, 7r+p, 
especially a separate measurement of the Dalitz contribution (e. g. with TAPS), could be 
used to verify the different realisations of the VMD. 
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TABLES 





a 


(MeV)- 1 


p + Be (1994) 
Pb + Au (1996) 


0.043 
0.035 


0.053 
0.023 



TABLE I. Parameters of the CERES momentum resolution (see text). 
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FIGURES 




M (GeV) 



Fig. 1 'Ideal' dilepton spectrum of p + Be collisions at 450 GeV. The 
UrQMD results are not filtered and the data points are for orientation 
only. 
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Fig. 2 Dilepton spectrum of p + Be collisions at 450 GeV. UrQMD 
calculation, corrected for acceptance and resolution, in comparison to 
data 
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Fig. 3 Dilepton spectrum of Pb + Au collisions at 158 A GeV. The '96 CERES 
data compared to filtered UrQMD calculations. 
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Pb (158 A GeV) Au 



2.1 <r;<2.65 
O ee > 35 mrad 
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Fig. 4 Effect of the mass dependence of the p meson dilepton decay width. Thin 
lines represent the mass distributions in the p channel while fat lines are the total 
spectra. (See text for more details.) 
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Fig. 5 Same as fig. |3| but with additional cuts on the pair trans- 
verse momentum qf +e of the dilepton. The low pt set is displayed 
in the upper part. 
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